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Abstract. We calculate for the nearest active galactic nucleus (AGN), Centaurus A, 
the flux of high energy cosmic rays and of accompanying secondary photons and 
neutrinos expected from hadronic interactions in the source. We use as two basic 
models for the generation of ultrahigh energy cosmic rays (UHECR) shock acceleration 
in the radio jet and acceleration in the regular electromagnetic field close to the core 
of the AGN. While scattering on photons dominates in scenarios with acceleration 
close to the core, scattering on gas becomes more important if acceleration takes place 
along the jet. Normalizing the UHECR flux from Centaurus A to the observations 
of the Auger experiment, the neutrino flux may be marginally observable in a Ikm^ 
neutrino telescope, if a steep UHECR flux dN/dE oc with a = 2.7 extends down 
to 10^^ eV. The associated photon flux is close to or exceeds the observational data 
of atmospheric Cherenkov and 7-ray telescopes for a >2. In particular, we find that 
already present data favour either a softer UHECR injection spectrum than a — 2.7 
for Centaurus A or a lower UHECR flux than expected from the normalization to the 
Auger observations. 



FAGS numbers: 98.70. Sa, 95.85.Ry, 95.85.Pw, 98.54.Cm 
1. Introduction 

Progress in cosmic ray (CR) physics has been hampered for long time by the 
deflection of charged cosmic rays in magnetic fields, preventing the identification of 
individual sources. Numerous searches for anisotropies and correlation studies have 
been performed, without reaching unanimous conclusions [1] . Using neutral messengers 
that should be produced as secondaries in proton-photon and proton-proton interactions 
close to the source for the identification of the sources has its own problems: First, 
secondary photons generated by hadronic CR interactions are difficult to disentangle 
from photons produced by synchrotron radiation or inverse Compton scattering of 
electrons. Moreover, high energy photons are strongly absorbed both in the source 
and propagating over extragalactic distances. By contrast, the extremely large mean 
free path of neutrinos together with the relatively poor angular resolution of neutrino 
telescopes (~ 1°) and the small expected event numbers makes the identification of 
extragalactic sources challenging using only the neutrino signal. Performing neutrino 
astronomy beyond the establishment of a diffuse neutrino background requires therefore 
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most likely additional input, either timing or angular information from high energy 
photon or CR experiments. 

The recently announced evidence [2] for a correlation of the arrival directions of 
ultrahigh energy (UHE) CRs observed by the Pierre Auger Observatory (PAO) with 
active galactic nuclei (AGN) may provide a first test case for successful "multi-messenger 
astronomy." In particular, Ref. [2] finds two events within the search bin of 3.1° around 
the nearest active galaxy, Centaurus A (Cen A or NGC 5128), that is located close to the 
supergalactic plane. At present this correlation has only 3 a CL. and other source types 
that follow the large-scale structure of matter would also result in an excess of events 
along the supergalactic plane. Independent evidence for the AGN source hypothesis is 
the characteristic bias of AGN with respect to the large-scale structure that is indeed 
reflected in the angular distribution of the observed UHECR arrival directions [3]. On 
the other hand, the correlation with AGN is not confirmed by the data from the HiRes 
experiment [4]. It is therefore timely to study the potential of high energy neutrino 
and photon observations for scrutinizing the correlation signal suggested by the Auger 
collaboration. 

The idea that neutrinos and photons are produced as secondaries in CR interactions 
close to the core of AGN has a long history [5]. In particular the diffuse neutrino and 
photon fluxes from all AGN have been studied in great detail [6-8]. Moreover, the 
expected neutrino and photon fluxes from Cen A were discussed recently in view of the 
PAO results in Refs. [9,10] and [11], respectively. The present work extends these later 
studies by calculating both neutrino and photon secondary fluxes and by modeling the 
source and particle interactions in more detail. Acceleration of protons to energies as 
high as 10^° eV in Cen A is used as an assumption that we discuss only as far as the 
source parameters are concerned. 

2. Source and acceleration models for Cen A 

2.1. Source parameters, the primary photon field and the gas column density 

Accretion A general review of the properties of Cen A, that is classified as PR I radio 
galaxy and as Seyfert 2 in the optical, is given in Ref. [12]. Its observed spectral 
energy distribution (SED) of electromagnetic radiation is discussed in detail by Ref. [13]. 
Although Cen A is the nearest active galaxy with its distance:}: oi D = 3.8 Mpc [14], 
there are a number of difficulties in deducing the photon distribution n^(x, p) that serves 
as target for hadronic interactions from observational data. First, the resolution of e.g. 
X-ray satellites like Chandra (0.5") or XMM-Newton (5") is not sufficient to resolve the 
core of the AGN. Second, part of the emitted radiation is heavily absorbed by the dust 
lane hiding the AGN core. Third, it is useful to distinguish between primary photons 
that serve as targets for hadronic interactions and secondary photons produced therein. 
Finally, the observed photon flux is at most energies probably dominated by photons 

X At this distance 1" = 18 pc. 
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produced in purely electromagnetic processes. 

We prefer therefore to model the primary photon field around the AGN core 
guided by the simplest possible theoretical model [15]. The thermal emission from 
a geometrically thin, optically thick Keplerian accretion disc is described by the 
temperature profile 

Within this simple model, the mass M of the central supermassive black hole (SMBH), 
its accretion rate M and the inner edge of the accretion disc Rq fix the main part 
of the primary electromagnetic radiation. Recent estimates for M vary between 
M = (0.5 - 2) X IO^Mq [16]. We will use M = 1 x lO^M© and thus the Schwarzschild 
radius is -R^ = 3 x 10^'^ cm. The angular momentum of the SMBH in Cen A is not 
known, so we use as smallest radius of the acceleration and of the emission region the 
radius of the last stable orbit for a Schwarzschild black hole, Rq = 3Rs ~ 1 x 10^^ cm. 

The accretion rate M and accretion efficiency r] = L/{Mc^) were fitted in Ref. [17] 
to Chandra and XMM-Newton observations assuming spherical (Bondi) accretion. The 
obtained accretion rate M = 6 x 10~^MQ/yr should be considered as a lower limit, 
because it does not account for the accretion of gas too cool to emit X-rays. In 
particular, it has been argued in Refs. [19, 20] that the heavy X-ray absorption and 
the high metallicity indicate accretion of cold gas. 

Integrating the surface brightness D = aT"^ over the finite accretion disc for the 
chosen values of M, M and Rq reproduces the characteristic blue bump [21]. Because of 
the relatively small accretion rate, the bump is shifted somewhat to lower frequencies 
compared to the standard case. The dust lane of Cen A prevents that such a bump can 
be seen by us directly in the SED of Cen A, but the observed HII line emission that is 
easiest explained by UV irradiation from the the central nucleus is indirect evidence for 
its existence [22]. 

Since the surface brightness drops fast with the radius, D oc r~^, most of the 
radiation is emitted close to the core, (3-15) cf. also Ref. [23]. To simplify our 
simulation, we consider therefore the following one-dimensional model for the source: We 
describe the accretion disc as a sphere of radius -Ri = 15Rs filled with an homogeneous, 
isotropic photon field radiating photons with the same spectrum n^{e) from each point 
on a "photosphere" with radius -Ri [24] . We model the energy-dependence of by 

{25£:cV £ev < 0.2 

e^Y 0.2 < £ev < 5 (2) 

0.l£^vexp(-eev/2) £eV > 5 

with ^ev = s/eV. The exponential cutoff is connected to the maximal temperature of 
the disc close to Rq. Finally, we assume that a hot corona produces an additional X-ray 
component, 

nx{e) = Kxs'y'^ , (3) 
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where the exponent —1.7 is chosen to agree with the X-ray observations from 
Refs. [17,20] and we use as high-energy cutoff e = 100 keV. 

Normalization The accretion rate M = 6 x lO~^M0/yr determines together with the 
accretion efficiencies rji the luminosity Li in the wave-length range i, Li = rji Mc^. For 
the X-ray range, we use Lx = 4.8 x lO^^erg/s in the 2-10 keV range according to the 
observations [17, 20] together with 

Li = ttRI c J de en^{e) = r/jLBondi = r/jMc^ , (4) 

giving Kx = 8 x 10^"*^/ (cm'^ eV) and rjx = 1.5%. The 7-ray luminosity = 5 x 10^^ erg/ s 
observed by COMPTEL and OSSE [25] corresponds to rj^ = 15%. Finally, we choose 
the normalization of the UV bump, Kyy in Eq. (2) as ?7uv = 10%. 

The used numerical values of the various normalization constants are summarized 
in Table 1. There we report also the resulting efficiencies f] = Li/LEdd relative to the 
Eddington luminosity, L-^dd = 1-5 x lO^^erg/s, of a 10*^ Mq black hole. The combined 
value of J2i Vi is in the range of advection-dominated accretion and thus the formation 
of a geometrically thin, optically thick accretion disc for Cen A is not guaranteed. We 
will nevertheless assume the existence of such a standard Shakura-Sunyaev accretion 
disk, keeping in mind however that the derived values for the photon density should be 
considered as upper limits. 

Table 1. The luminosities, efficiencies and normalization constants in different wave- 
length ranges. 



band 


Li/erg/s 


V 


V 


Ki/ {cm^ eV) 


7-ray 


5.0 X 10^2 


15% 


3.3 X 10-^ 




X-ray 


4.8 X 10^^ 


1.5% 


3.3 X 10"^ 


8 X 10" 


UV 


3.6 X 10^2 


10% 


2.4 X 10-^ 


1 X 10^3 



For large distances, r ^ Ri, the photon field emitted from the photosphere of radius 
Ri scales as n-y(r) oc (i?i/r)^ and has approximately a distribution of momentum vectors 
with cos'd = p- r/(|p||r|)>l — (i?i/r)^. Thus we have to distinguish two cases: Either 
photons and/or cosmic rays are distributed isotropically or both are non-isotropic. The 
latter case is realized when charged particles propagate along the (regular) field lines 
and primary photons stream nearly radial outwards for r ^ Above threshold 

E ~ 10^^ eV, the resulting interaction depth is Tp^ ~ 3. 

The PAO data used in the analysis [2] correspond to an exposure S = 9000 km^ yr sr 
with maximal zenith angle "(^max = 60°. The correlation signal with AGN was maximized 
for the threshold energy i?th = 6 x 10^^ eV and the angular bin size 3.1°. For these values, 
two events were found in the angular bin around Cen A. Centaurus A is close enough 
to the Earth that energy losses of CRs with energy close to Eth can be neglected, cf. 
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e.g. Fig. 1 of Ref. [26]. Assuming that both events indeed originate from this AGN, the 
integral CR flux above E^ii on Earth from Cen A is 

F(> = 1^ = 2 X 10-3 . (5) 

R ^ km^ yr 

Here, ^ 9 sr is the field-of-view of the PAO and R = rj{6s)/ {f]{^)) ~ 0.95 is the ratio 
of the exposure rjlSg) at the declination 6s = —43.0° of Cen A and the average exposure 

iv)- 

Gas column density and proton-proton interactions Reference [20] fitted Suzaku 
observations of Cen A with several absorption models. The primary X-ray component 
was found to be absorbed by the column density X = Ix lO^^/cm^, while a less brighter 
component is more heavily absorbed with X = 7 x lO^^/cm^. These results can be 
interpreted either as indication for a clumpy structure of the gas or for two different 
X-ray sources, e.g. from accretion and from jet emission. The resolution of these 
observations is however rather limited, and the value for the column density is therefore 
biased by the dense torus. A more representative determination of the mean density 
of the gas around the core of Cen A was possible with XMM-Newton and Chandra 
observations. Ref. [27] determined the density of the interstellar medium around the 
core of Cen A as uh = no[l + (r/ro)^]"*^'® with uq = 0.04/cm3 and tq = 0.5 kpc. 
Finally, Ref. [28] found as average hydrogen density X = 1.5 x 10^^/cm^ along the 
radio jet starting from a projected distance of 0.3 kpc up to 2.5 kpc. With d = 0.4 kpc 
as diameter of the jet and assuming that the jet axis is close to perpendicular to the 
line-of-sight [29], an average density uh ~ Xu/d ~ 1.7/cm'^ follows. 

2.2. CR acceleration, propagation and interactions in the source 

X-ray observations of Cen A allow one to trace the acceleration sites of electrons 
along the radio jet, because their synchrotron loss length is short compared to the 
extension of the jet. The strong dependence of synchrotron emission on the mass of the 
radiating particle implies that the SED is, apart from the highest energies, dominated 
by electromagnetic interactions of electrons and that therefore acceleration site and 
mechanism for electrons and protons may differ. Although various proposal for the 
acceleration of protons to UHE exist, we restrict ourselves to two basic models, namely 
shock acceleration in the radio jet and acceleration in the regular electromagnetic field 
close to the core of the AGN. These two models are characterized by a rather different 
set of parameters and the secondary fluxes in several other models may be obtained 
by "appropriate interpolation." We neglect relativistic effects, because of the moderate 
Lorentz factors observed in Cen A and the large angle between the jet axis and the 
line-of-sight. 

Acceleration in regular fields near the core Acceleration close to the core could proceed 
either via shock acceleration in accretion shocks [30] or via acceleration in regular 
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fields. The former case is disfavoured since the acceleration rate is smaller than the 
rate of photonpion energy losses [31]. Therefore we consider here only acceleration 
in regular fields close to the core. Among the possibilities discussed in the literature 
are electromagnetic winds from the BH/disk magnetosphere [32, 33], and unipolar 
induction around a Kerr BH [34, 35]. The regular magnetic field close to the core 
consists of a toroidal component in the accretion disc and a poloidal field component 
Bp = BQ{Ro/r)^ . For field strengths Bq ~kG acceleration to energies around (10^'' eV 
is possible. Since the curvature radius Rc of the field-lines is typically large, curvature 
radiation of protons is not very effective. Moreover, protons move mainly parallel to the 
field-lines and thus also synchrotron losses are suppressed. Because of our simplified one- 
dimensional geometry, we assume that the combined effect of synchrotron and curvature 
radiation is such that is does not affect protons but acts as main energy loss for electrons. 
The other remaining free parameter in this model, the injection point of the CRs, is 
fixed to r = GRg. 

The required magnetic field strength for acceleration to 10^'' eV is an order of 
magnitude higher than one would expect from equipartition, B^/Sn ~ L/AttR^c. The 
currently low accretion rate and small luminosity of Cen A disfavor therefore acceleration 
close to the core as acceleration mechanism, if Cen A was not in the past in a state of 
higher activity. 

Acceleration in the radio jet Rachen and Biermann suggested the hot spots in FR- 
II radio galaxies as sites of UHECR acceleration [36]. These spots are formed as 
termination shock of the supersonic jets in the intergalactic medium and are especially 
prominent in FR II jets. They are instead dim or absent in weak FR I sources, most likely 
because strong turbulent dissipation in the propagation phase reduces the momentum 
finally released at the termination shock. Nevertheless, e.g. Ref. [37] argued that at the 
"hot spot" in Cen A acceleration of protons to UHE occurs. The projected size of the 
hot spot is Rhs ~ 1.7 kpc, containing magnetic fields with strengths B ^ 0.5 mC. 
Alternatively, proton could be accelerated along the whole extension of the radio 
jet, extending a projected distance of 6 kpc from the nucleus. In this case, protons 
have to diffuse through the X-ray photons emitted by accelerated electrons and, more 
importantly, through the hydrogen gas observed in the jet. 

Thus we shall use acceleration in the jet as our second basic scenario. More precisely, 
we consider as acceleration region for protons a cylinder of length / = 4 kpc and diameter 
d = 0.3 kpc, similar to the emission volume of X-rays observed by Chandra. Diffusion in 
the turbulent magnetic fields will increase the interaction depth. We choose as field 
strength B = 0.2 mC, set the coherence length conservatively equal to the length 
Ic = Ikpc and use a = 1/2 for the energy dependence of the diffusion coefficient, 
D[E) oc corresponding to a Kraichnan fluctuation spectrum. The choice a = 1/2 
is intermediate between the often assumed Bohm and Kolmogoroff diffusion as well as 
close to the numerical results of Ref. [38] . 
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Relative importance of photon and proton targets In the case of acceleration close 
to the core, either in accretion shocks or regular electromagnetic fields, UV photons 
are the most important scattering targets and the interaction depth for photo-hadron 
interactions can reach Tp^ ~ few. In all other will be reduced at least by a 

factor ~ lORs/l where I is the typical dimension of the acceleration region. Diffusion 
increases the effective size of the source, but this effect becomes noticeable mainly at 
energies below the threshold for photon-hadron interactions. As a rule of thumb, photo- 
hadron interactions are therefore only important when the acceleration takes place close 
to the core. 

The importance of proton-proton interactions depends more strongly on the 
concrete model assumptions - as the scatter in the observed column density X of gas 
shows. Moreover, the accretion rate may be time-dependent and matter close to the core 
is concentrated in clumps. Thus both the CR luminosity Lqr and the interaction depth 
Tpp can vary significantly depending on the model assumptions. While the importance 
of gas as target close to the core is uncertain, the Chandra observations indicate that 
pp scattering is the main source of CR interactions in the jet. 



Energy spectrum and interactions We consider three cases for the generation spectrum 
dN/dE oc E~'^ of UHECRs: i) A power-law with a ^ 2 as conventionally predicted 
from first-order Fermi acceleration with non-relativistic shocks, ii) A broken power- 
law with break energies Ei, = 10^'^ eV and Ei, = 10^^ eV, respectively, and a = 2.7 for 
E > Ef, as suggested by the dip-interpretation of the experimental data [39]. Since the 
steepening of the observed diffuse spectrum may result from a distribution of maximal 
energies dn/d£'max [40], present UHECR observations do not distinguish between these 
two possibilities, even if one assumes that extragalactic sources dominated the flux down 
to 10^'' eV. Note also that a = 2.7 is consistent with the observed radio spectrum of 
most AGN. In order to weaken the energy demands of Cen A, we extend the a = 2.7 
spectrum with a = 2 below E^. Hi) A flat spectrum with a = 1.2 and most energy 
concentrated at -Emax, as expected for a "linear accelerator." In all three cases we use 
E = 10^° eV 

Hadronic interactions are simulated with an extension of the Monte Carlo code 
described in Ref. [41], including the possibility of a non-thermal and anisotropic photon 
fleld n^(x, p) and of proton-proton interactions. Electromagnetic processes§ like pair 
production and inverse Compton scattering of electrons and photons both inside the 
source and on IR, cosmic microwave background (CMB) and radio photons are also 
taken into account. Finally, diffusion is simulated as described in Ref. [41]. 



3. Cosmic ray, photon and neutrino fluxes 

The left and right panels of Fig. 1 display the particle fluxes F{E) predicted from Cen A 
as function of the energy E in the case of acceleration close to the core and in the jet, 

§ A description of our treatment of electromagnetic cascades and results will be presented elsewhere. 
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Figure 1. Particle fluxes F{E) from Cen A normalized to the PAO data as function 
of the energy E: top panels broken power-law with a = 2 and 2.7 and break energy 
Eb = 10^^ eV, middle a — 2, and bottom panels a — 1.2; left panels acceleration close 
to the core, right panels in the jet; initial protons (solid black), final protons (dashed 
black), photons (blue), and sum of all neutrinos (red). 
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respectively. Additionally to the injected proton flux (black solid line), we show the 
flux of protons (black dashed) arriving on Earth. The flnal proton flux is reduced by 
interactions by a factor 2 above the threshold energy ~ 10^^ eV for photon-proton 
interactions (left), while diffusion in the jet increases the interaction depth for lower 
energies (right panels), resulting in the effective production of secondaries. 

The photon flux at Earth after cascading in the source and on the extragalactic 
background light (EBL) consists of two contributions: i) Photons produced in the 
source, which survive their subsequent travel to the observer without interactions with 
the EBL; ii) photons produced during the cascade on the EBL. While the former are 
always observed as a pointlike contribution from a pointlike source, the latter are 
characterized by a flnite angular spread, due to the deflections of cascade electrons 
in the extragalactic magnetic fleld (EGMF). Depending on the strength of the latter 
and on the experimental resolution, the second contribution may be observed as a halo 
around the source [45], while contributing only partially to the unresolved pointlike flux. 
Because of the proximity of Cen A, deflections in the EGMF have a — compared to the 
overall uncertainties — negligible influence on the calculated photon spectra. 

The photon flux at Earth after cascading in the source and on the ELB is shown 
as a solid blue line together with the combined flt to EGRET, OSSE and COMPTEL 
observations of Cen A from Ref. [25], the H.E.S.S. limit from Ref. [46] for a = 2 and 
a = 3, and the FERMI sensitivity [47] for a 5cr detection of a point source in 1 yr. The 
predicted spectrum shows the typical suppression above the pair production threshold 
on the CMB at ~ 200 TeV. Since the CR spectra are normalized to the integral 
UHECR flux above i^th = 5.6 x 10^^ eV, steeper spectra result in larger secondary fluxes 
at low energies. In particular, the photon flux overshoots the H.E.S.S. limit or the 
CGRO observations in the case of a broken power-law injection spectrum. Thus already 
present observational data favour either softer UHECR spectra than dN/dE oc 
for Cen A or a lower UHECR flux than expected from the normalization to the Auger 
observations. Note however that our normalization relies on only two events and has 
therefore a large statistical uncertainty. On the other hand, the determination of the 
energy scale of UHECR experiments is notoriously difficult and it has been argued that 
the PAO energy scale should be shifted up to obtain agreement with the spectral shape 
predicted by e^e~ pair production [39]. 

Comparing the neutrino flux at 10^^ eV for the different models allows one to 
judge how strong the predicted neutrino event number depends on the slope of the 
CR spectrum. Going from a broken power-law with a = 2.7 at UHE to a = L2 reduces 
between three and four orders of magnitude the neutrino flux at 100 TeV. Calculating 
the expected event number in a neutrino telescope requires a deflnite choice of the 
experiment. Centaurus A is from the location of Icecube only visible from above, and 
thus the background of atmospheric muons allows only the use of contained events 
that carry essentially no directional information [6'& ~ 30°). For the calculation 
of the number of contained events expected in Icecube we use as effective volume 
V = Ikm^, as threshold -Emin = 100 TeV, assume 100% efficiency above -Emm and use 
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the CTEQ5 neutrino-nucleon cross sections [42]. By contrast, a neutrino telescope in 
the Mediterranean could make use of the muon signal and the directional information. 
In this case the rate R of muon events can be approximated by 

R = A dES{E)F,^{E)P{E,E^^) (6) 

with A = 1 km^, the probability -Emm) = Na_{R{E, -E'mm))o"i/Ar that a muon reaches 
with E > -Emin the detector and the angular averaged shadowing factor S{Ey) accounting 
for attenuation of the neutrino flux in the Earth [43]. The resulting event numbers both 
for cascade and shower event number per year observation time are summarized in 
Table 2 together with the input CR luminosity. Note that the cutoff in the neutrino 
spectra below 100 GeV that is visible in the upper, right panels of Fig. 1 is artificial, 
since we neglect neutrinos with lower energies in our simulation. 

In summary, we predict that for all cases where there is a marginal chance to detect 
a neutrino signal from Cen A, atmospheric Cherenkov telescope and/or 7-ray satellites 
like FERMI should detect previously a photon signal from Cen A. Main reason for this 
result — that is in contradiction to earlier expectations, cf. e.g. Ref. [44] — is the cascading 
of photons in the anisotropic photon field close to the source. 

Finally we remark that we assumed for the calculation of all fluxes an isotropic 
emission. This is well justified in the case of stochastic acceleration in the jet, because 
of the small gamma factors of the observed matter flows in the jet and the large angle 
between the jet and the line-of-sight. In the case of acceleration close to the core, 
protons move along the field lines and thus the emission is rather anisotropic. Since 
we used the UHECR flux towards our line-of-sight as normalization, the total UHECR 
luminosity of Cen A is therefore larger than the one estimated in Table 1. On the other 
hand, the predicted number of neutrino events is not affected by the anisotropy, since 
the ratio of neutrino and UHECR fluxes is roughly independent from the considered 
direction. By contrast, TeV photons as final products of electromagnetic cascades are 
more effectively isotropized than UHECRs. This effect will slightly increase the TeV 
photon flux compared to the isotropic case assumed in Fig. 1. 



Table 2. The injected CR luminosity LcR and the number of neutrino events expected 
per year observation time for different energy slopes and acceleration scenarios. 



a or Eb/eV 


jet 


core 


1.2 


2.0 




10^^ 


1.2 


2.0 




10^^ 


LcR/lO^Oerg/s 


0.35 


0.97 


5.2 


5.2 


1.1 


2.2 


11 


10 


contained ^ u/ji 


8 X 10-^ 


0.02 


0.4 


2.0 


7 X 10-^ 


0.01 


0.3 


0.9 




4 X 10"^ 


7 X 10-3 


0.2 


0.7 


3 X 10"^ 


7 X 10-3 


0.1 


0.5 



The expected neutrino and photon fluxes from Cen A were discussed recently in 
view of the PAO results in a rather qualitative way in Refs. [9, 10] and [11], respectively. 
Our case of acceleration close to the core and a broken power-law corresponds roughly 
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to the set-up of Ref. [9], but the results differ by the large factor 30. One possible 
explanation of this difference is that their neutrino flux drops faster than we found. 
Reference [10] assumed pp interactions on gas close to the core with Tpp ~ few as 
production mechanism. The obtained event numbers agree well taking into account the 
differences in the used assumptions. Finally, Ref. [11] discussed photons from hadronic 
interactions in Cen A. This work did not include the effect of electromagnetic cascading 
and the conclusions are therefore difficult to compare. 

4. Conclusions 

We have calculated the flux of high energy cosmic rays and of accompanying secondary 
photons and neutrinos expected from Cen A. We modeled the distribution of target 
photons and gas guided by the simplest theoretical model for the accretion disc and by 
observational data, respectively. The production of secondaries and the electromagnetic 
cascading of electrons and photons was simulated with a Monte Carlo procedure. 

In contrast to previous works, we showed that scattering on gas becomes important 
if acceleration takes place along the jet. Moreover we found that a source that has 
an interaction depth Tp^ > 1 can be observed in the (1-100) TeV range by atmospheric 
Cherenkov telescopes. Additionally to these more technical results, we have shown 
that a combination of the old CRGO observations and the limits from atmospheric 
Cherenkov telescopes can be used to constrain currently favoured UHECR models. In 
particular, we found that these data favour either a softer UHECR injection spectrum 
than dN/(\E oc i?"^'^ for Cen A or a lower UHECR flux than expected from the 
normalization to the Auger observations. 

One should remind however the main underlying uncertainties interpreting our 
results: The normalization of all fluxes is based on the assumption that two UHECR 
protons in the PAO data originate from Cen A. Note that several authors have argued 
that a larger number of events originates in Cen A [48]. Apart from the purely Poisson 
error, the normalization may be influenced by deflections in (extra-) galactic magnetic 
fields, the uncertainty in the energy scale of PAO, and a possible admixture of heavy 
nuclei. Deflections of UHECRs result in their delayed arrival with respect to photons and 
neutrinos, introducing an additional source of uncertainty in their relative normalization. 
Moreover, our model parameters (M, M, r/uv,- • • ) and even as basic parameters as the 
distance to Cen A have sizeable uncertainties. Last but not least we have made several 
simplifying assumptions like the use of an one-dimensional geometry and the omission 
of the acceleration process: We only postulated that acceleration to 10^° eV is possible 
in the environment of Cen A, without demonstrating it for a concrete model. Despite 
these drawbacks, it is remarkable that 7-ray and TeV observations of Cen A allow one 
already now to constrain currently favoured UHECR models. The potential of neutrino 
telescopes to observe Cen A depends strongly on the steepness of the UHECR generation 
spectrum. A neutrino telescope on the northern hemisphere would be very useful for 
this task. 
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